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The degree of electron correlations remains a central issue in the iron-based superconductors. The parent 
iron pnictides are antiferromagnetic, and their bad-metal behavior has been interpreted in terms of proximity 
to a Mott transition. We study such a transition in multi-orbital models on modulated lattices containing an 
ordered pattern of iron vacancies, using a slave-rotor method. We show that the ordered vacancies lead to a 
band-narrowing, which pushes the system to the Mott insulator side. This effect is proposed to underlie the 
insulating behavior observed in the parent compounds of the newly discovered (K,Tl)yFexSe2 superconductors. 



Introduction. Superconductivity in the layered iron pnic- 
tides and chalcogenides occurs near antiferromagnetically- 
ordered parent compounds In their paramagnetic 

phase, these parent materials have a large electrical resistivity 
corresponding to an in-plane carrier mean-free-path on the or- 
der of the inverse Fermi wavevector They also show a strong 
reduction of the Drude weight 101, and temperature-induced 
spectral weight transfer that extends to high energies (on the 
eV order) i5i(-|7l]. Such bad-metal behavior is characteristic of 
metallic systems in proximity to a Mott transition Isl- lioll 

Recently, superconductivity has been discovered in a new 
family of iron-based compounds KyFe2Se2 ifTlll and related 
(K,Tl)FexSe2 lll2ll . In these compounds the maximal super- 
conducting transition temperature is comparable to that of 
the 122 iron pnictides. Similarly to the pnictides, the super- 
conductivity occurs close to an antiferromagnetically ordered 
state 1121]. At the same time, these materials are unique in sev- 
eral aspects. Both the angle-resolved photoemission (ARPES) 
experiments lO-lH] and LDA calculations fl6ll show that 
the Fermi surface has only electron pockets. The absence 
of hole Fermi pockets is unique among the iron based su- 
perconductors, raising hope for major new and general in- 
sights to be gained from studying these materials. Equally 
important, the Fe vacancies may form ordered patterns when 
the Fe content x < 1.6 as suggested by various experi- 



ments rtl2lll7Lll8[l . Furthermore, there are parent compounds 
which are insulating 1 12 , 3 ■ The control parameter that tunes 
the (K,Tl)FexSe2 system from superconducting to insulating 
is the Fe composition x, and x ~ 1.5 is the primary candi- 
date composition for a parent compound. There is evidence 
1 1211171] that in (K,Tl)Fei.5Se2 the Fe vacancies form regular 
patterns possibly as illustrated in either Fig. [Ttb) or Fig. [IJc). 
The in-plane electrical resistivity is about two orders of mag- 
nitude larger than that of the parent iron pnictides at room 
temperature, and it further increases exponentially as temper- 
ature is lowered. The insulating behavior is also manifested 
in the optical conductivity 1I20I1 . which is strongly suppressed 
below about 0.7 eV. Because of experimental indications that 
the (K,T1) content is also variable, we will in the following 
refer to these systems as (K,Tl)yFexSe2. 

In this Letter, we propose that the parent (K,Tl)yFei 5Se2 
is a Mott insulator arising from a correlation effect that is en- 
hanced by the Fe vacancies. We describe the ordered Fe va- 



cancies in terms of a modulated lattice, and introduce a two- 
orbital model with two electrons per Fe site to capture their 
electronic structure. We use a slave-rotor method to show that 
a Mott transition exists in this model even though there are 
an even number of electrons per site (and per unit cell). We 
find that the interaction strength for the Mott transition largely 
tracks the electronic bandwidth. In other words, ordered Fe 
vacancies enhance the tendency towards Mott localization as a 
result of a kinetic-energy reduction. Such Fe vacancies, there- 
fore, have a similar effect as a lattice expansion, which we 
have previously discussed as responsible for the Mott insulat- 
ing behavior in La203Fe2Se2 12111 . Our considerations of the 
interaction effects are realistic, given that the ab initio calcu- 
lations using density-functional theory l22ll show that the 3d 
bands of TlFei.5Se2 are narrower than those of TIFe2Se2. 

Modulated lattices and kinetic-energy reduction. We will 
consider a square lattice (LI, Fig.[TJa)), a modulated square 
lattice consisting of 2 x 2 plaquettes each having its center 
removed (L2, Fig. [lib)), and another one coiTesponding to a 
triangular lattice of such 2x2 plaquettes (L3, Figfljc)). 

The (K,Tl)yFexSe2 system involves all five 3d orbitals. 
The ARPES experiments Il3l4l5ll show electron pockets near 
M point and suggest very weak electron-like pockets near F 
point. The absence of hole pocket near F point is largely con- 
sistent with the ab initio electronic bandstructure calculations 
usin g lo cal-density-approximation (LDA) for (K,Tl)Fe2Se2 



1 16l I23l - l26ll . This is in contrast to the case of iron pnic- 
tides, and is easier to model using a two-orbital tight-binding 
parametrization. CoiTespondingly, we consider a two-orbital 
model with the degenerate xz and yz orbitals (labeled as or- 
bitals 1 and 2) and n = 2. Inspired by the considerations in 



the pnictides case 127112811 . we introduce a set of tight-binding 
parameters. The parameters are listed in Table U and their 
meanings can be inferred from the dispersion functions speci- 
fied in Eq. (|2]i. We first fit the LDA bandstructure obtained on 
TlFe2Se2 to this two-orbital model and then adjust the tight- 
binding model parameters so that the Fermi surface still has 
only electron pockets at n = 2. We notice that the Fermi 
surface size is larger than in the bandstructure calculations, 
but this suffices for our qualitative considerations of the effect 
of lattice depletion on the Mott transition. What is important 
is that, for our parameters, the Fermi surface comprises only 
electron pockets near the X points of the 1-Fe per cell Bril- 
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louin Zone. The bandwidth narrowing for L2 and L3 lattices 
compared to the LI lattice is shown in Fig. [Tie)- 
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-0.038 
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TABLE I. Hopping parameters of the two-orbital model. 



Two-orbital model and the slave-rotor method. We are now 
in position to specify our model, 
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FIG. 2. (Color online) Quasiparticle weight for the one-orbital model 
in the unmodulated and modulated lattices plotted as a function of 
U/t (a) and U/D (b). 



where Ciaa annihilates an electron in orbital a and spin a on 
site i of the Fe lattice. The first term in Eq. ([T]i describes the 
electron hopping, with orbital dependent hopping amplitudes 
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FIG. 1. (Color online) a) Square lattice, LI. Different clusters are 
used and will be labeled as Nc = 1,2, and 4, respectively; b) Modu- 
lated 2x2 square lattice, L2. The enclosed sites 1,2,3 form the basis 
of the unit cell; c) Another modulated 2x2 lattice, L3 (correspond- 
ing to a 4 X 2 superstructure in the FeSe plane). Also shown are the 
bare density-of-states (DOS) for the one-orbital model (t = t' = 1) 
(d) and the two-orbital model (e). 



ttf = ti,t2, ti2 yielding T,kaf3a 4a4aaCkaa, whcre 

e]} = —2ti COS kx — 2t2 cos ky — At^ cos kx cos ky 

—Atg cos 2kx cos 2ky, 
ef^ = —2t2 cos kx — 2ti cos ky — At^ cos k^ cos ky (2) 



'Atg cos 2kx cos 2ky, 



VI _ ,21 



-4t4 sin kx sin ky — 4ti2 sin 2kx sin 2ky 



Small hoppings tg and t\2 between the 5th-nearest neighbors 
are included to reproduce the two electron pockets at n — 2. 
The second term in Eq. (|T) is an on-site Coulomb repulsion. 
We focus on the effect of lattice depletion on the Mott transi- 
tion and will not consider Hund's coupling and pair-hopping 
terms for simplicity. All local interactions are expected to 
have effects similar to I] . In particular, the Hund's coupling 
will reduce the critical U of the Mott transition; its effects are 
readily studied within a slave-spin method ll29[ Isoll . and the 
results will be reported elsewhere. 

We study the model using the cluster slave-rotor mean- 
field (CSRMF) method fHiil. We introduce an 0(2) ro- 
tor variable Qi and a spinon Jiaa on each site, and write 
Ciao — /iQcre"'^'. Here, e~*^' lowers the rotor angular mo- 
mentum Li, which corresponds to the charge quantum num- 
ber The unphysical states are eliminated by enforcing the 
constraint = J2aMLafiacr - 1/2) in the enlarged rotor 
and spinon Hilbert space. By rewriting Eq. ([T]i using rotor 
and spinon operators and decoupling the rotor and spinon op- 
erators at the mean-field level, we obtain the following two 
effective Hamiltonians: 



J 



(3) 



where C, 



) ) fl is the rotor correlation function that 



renormalizes the quasiparticle hopping parameters in the pres- 
ence of interaction, xl^ = {fiaafjl3<^)f^ M is the chemical 
potential, and A is a Langrange multiplier to impose the con- 
straint. To solve these two Hamiltonians, which still contain 
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FIG. 3. (Color online) Quasiparticle weight Z for the two-orbital 
model vs. I] jt (a) and U /D (b). Also shown are the site-selective 
Zi and Zo, and the bond correlations for lattices L2 (c) and L3 (d). 



interactions among rotors, we further apply a cluster mean- 
field approximation. We exactly diagonalize the rotor Hamil- 
tonian on a finite cluster, and treat the influence of the sites 
outside the cluster as a mean field. We decouple e*'^'^^^' 
into e^^^(j>j if i belongs to the cluster but j is outside. Here 
= (e~*^') is the local mean-field parameter Consequently, 
Cij is factorized as w 4'i4'j if either i or j is outside 
the cluster In practice, Eqs. Q and (|4|i are self-consistently 
solved by iteratively determining the mean-field parameters 
(/){, Cij, and xtf ■ Mott transition is signaled by a vanish- 
ing quasiparticle spectral weight Z = \4>i\'^. 

To gain intuition on the role of the lattice modulation, we 
will also study a one-orbital model with nearest-neighbor hop- 
ping, t, and next-nearest-neighbor hopping, t' . Fig.fljd) illus- 
trates the reduction of the bandwidth for L2 and L3 lattices 
from that of the LI lattice, for the case of t' = t. A non-zero 
t' is chosen for two reasons. It avoids a perfect nesting in the 
case of n = 1, which we study below. It also avoids a flat band 
in the case of the L2 lattice: when t' = 0, for the L2 lattice, 
there are two dispersive bands with a combined bandwidth of 
4:^/2t, and a flat band in the middle. 

Results for the one-orbital model. We start from the one- 
orbital case. Because the L2 and L3 lattices involve a 2 x 2 
square plaquette as the unit cell, we will carry out our calcu- 
lations for the lattice LI with Nc = 4. The slave-rotor mean- 
field theory treats the rotor kinetic energy for intra-cluster 
bonds exactly by diagonalizing the rotor Hamiltonian on the 
cluster. Hence working with Nc = ^ gives a better description 
of the Mott transition than using the single site approximation. 
Fig. 12 a) shows the renormalized quasiparticle weight, Z, as 
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FIG. 4. (Color online) Sketch of the renormalized hoppings for the 
L2 lattice at [/ < Ud (a) and at Uci < U < Uc2 (b). Also shown 
are the bandstructure of L2 for U — < Ud (c) and for Uci < U = 
19ti < Uc2 (d). The bands are shown along the high-symmetry 
directions in the Brillouin Zone associated with the 2 x 2-Fe unit 
cell, illustrated in the inset in (d). 



a function of U /t. The Mott transition occurs at Uc, where Z 
first goes to zero as U is increased. For the LI lattice, increas- 
ing the cluster size from iVc = 1 to iVc = 4 leads to a succes- 
sive reduction of Uc'- Uc ~ 8t for Nc = 1, and Uc ~ 7.2t for 

A^c = 4. 

On the modulated lattices, we also find a Mott transition. It 
is seen that Uc ~ 5.4i for the L2 lattice and Uc ~ 5.8t for L3 
lattice; both are smaller than that of the LI lattice. 

Fig.|2|b) plots the same result, but with U now normalized 
against the full bandwidth D. It is seen that Uc/D is compa- 
rable for all three cases. This clearly illustrates that the reduc- 
tion of Uc for the modulated L2 and L3 lattices arises from the 
band-narrowing effect. 

Results for the two-orbital model. We now turn to the 
more realistic two-orbital case. The renormalized quasipar- 
ticle weight as a function of U /ti is shown in Fig. [3ja). It is 
again seen that the values of Uc for both the L2 and L3 lattices 
are smaller than that of the LI lattice with Nc = 4. 

One difference from the toy 1 -orbital model is that the 
hopping parameters in the two-orbital model are highly 
anisotropic (e.g., |t3/i2| ~ 3), which makes the local envi- 
ronment possibly different from site to site on the modulated 
lattices. To fully address the influence of the inhomogeneity, 
we study the quasiparticle weight associated with each site 
in the cluster Zi — |(e*''')p. For either LI or L3, we ob- 
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tain a single Mott transition as in the one-orbital case. For 
the L2 lattice, we find two transitions. They are identified by 
the vanishing of Z3 first at Ud/ti « 17, and subsequently 
the vanishing of Zi (and, equivalently, Z2) at a higher value 
Uc2/ti « 20. 

To understand this, we have in the same figure plotted the 
bond correlators C12 and C13. Between Uci and Uc2, C13 
vanishes but C12 remains finite. This makes site 3 to be un- 
connected to the rest of the lattice {cf. Fig. Etb)). We see 
these explicitly in a plot of the renormalized bandstructure in 
Fig. lUd): associated with the isolated 3 sites is a flat band 
lying exactly on the Fermi level for Ud < U < Uc2- By con- 
trast, for U < Uci, all sites are connected by hopping terms 
(Fig.Uta)) and there exists no flat band (Fig.HJc)). 

For the L3 lattice, the geometry prevents the separation of 
any site from the rest bulk unless all the effective hopping 
parameters are zero. As a result, there will be only one transi- 
tion. This is clearly seen in comparing Zi and Z3 in Fig.[3jd). 

The quasiparticle weight as a function of U/D is shown in 
Fig. Ob). It is again seen that Uc/D is comparable for all 
three cases. As in the one-orbital case, this illustrates that the 
reduction of Uc for the modulated lattices originates from the 
band-narrowing effect. 

Implications for (K,Tl)yFe^Se2. Our results imply that 
the critical interaction strength for the Mott transition will be 
smaller in (K,Tl)yFexSe2 than in iron arsenides and 1 1 iron 
chalcogenides. This provides the basis for a Mott-insulating 
state in the parent (K,Tl)yFexSe2, even when one assumes the 
same strength of Coulomb interaction across all families of 
iron-based superconductors.. 

The Mott-insulating nature of the parent (K,Tl)yFexSe2 is 
supported by experiments. As already mentioned, the materi- 
als for both X = 1.5 and x = 1.64 have a large electrical re- 
sistivity with an insulating temperature dependence lll2lll9[l .. 
Furthermore, the insulating behavior in the electrical resistiv- 
ity is already observed in the paramagnetic phase. Relatedly, 
the optical conductivity is not only strongly suppressed be- 
low about 0.7 eV, but also small in magnitude. For reference, 
the value of the optical conductivity is comparable to that of 
the insulating YBa2Cu306+i^ with a slight off-stoichiometry 
X = 0.2 [330 • Finally, magnetic order is known to exist in 
TlFexSe2 at x close to 1.5 l34ll . Taken together, these exper- 
iments suggest that the insulating state is of the Mott type. 

We note that in compounds with Fe content close to x_ 
1.6, the ordered vacancies have a different pattern 1112 , 



expected on the ground of general considerations given here, 
and can also be seen in the LDA results 1220 . 

To summarize, we have used a two-orbital model in 1/4- 
depleted lattices to demonstrate that ordered vacancies en- 
hance the tendency for Mott transition, and that this enhance- 
ment originates from a vacancy-induced kinetic-energy reduc- 
tion. Our qualitative conclusion is expected to apply to the 
more realistic five-orbital model. Based on our calculations, 
we propose that the insulating parent of the (K,Tl)yFexSe2 
superconductors is a Mott insulator at ambient pressure. 
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sions. This work was supported by NSF Grant No. DMR- 
1006985 and the Robert A. Welch Foundation Grant No. C- 
141 1 (R. Y. and Q.S.), and U.S. DOE at LANL under Contract 
No. DE-AC52-06NA25396 (J.-X.Z.). 



However, because we have shown that the Mott localization 
is a result of vacancy-ordering induced band narrowing, our 
argument will also apply to these systems. Band narrowing is 
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